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Abstract 

At carrier densities above the Mott density Coulomb screening destroys the exciton resonance. 
This, together with band-gap renormalization and band filling, severely affects the optical spectra. 
We have experimentally studied these effects by ultrafast pump-probe reflectivity measurements on 
a ZnO single crystal at various wavelengths around the exciton resonance and in a broad carrier- 
density range. Theoretically we determined the Mott density in ZnO to be 1.5 x 10 24 m -3 at 300 
K. Taking a field-theoretical approach, we derived and solved the Bethe-Salpeter ladder equation 
and we computed the density-dependent reflectivity and absorption spectra. A carrier dynamics 
model has been developed, containing three-photon absorption, carrier cooling, and carrier trapping 
near the surface. The agreement between the theoretical reflectivity based on our model and the 
experimental data is excellent. 

PACS numbers: 78.47.jg, 78.20.Bh, 79.20.Ws, 71.35.-y 
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I. INTRODUCTION 



Because of its wide direct band gap of 3.37 eV, ZnO has many possible applications 
for optoelectronic devices, including solar cells and light-emitting diodes. ZnO nanowires 
are used as waveguides and UV lasers,- - - photodetectors,- and optical switches.- For such 
applications, it is important to know and understand the optical spectra of ZnO at high 
carrier densities, as well as the carrier dynamics. 

The optical spectra at high carrier densities are strongly influenced by screening, band-gap 
renormalization (BGR), and band filling. At densities higher than the so-called Mott density 
hm, screening of the Coulomb interaction destroys the exciton resonance. Here we present a 
concise theoretical and experimental study of these phenomena in ZnO, covering the exciton 
regime, the electron-hole plasma regime, and the crossover between them. Analysis of the 
pump-probe reflectivity experiment described here also reveals the ultrafast carrier dynamics 
near the crystal surface. 

Pump-probe reflectivity experiments on ZnO by other groups 7-9 have shown large reflec- 
tivity changes at high carrier densities. Despite these striking observations, a quantitative 
picture of the carrier dynamics and the reflectivity spectrum at high densities does not 
exist. This has several causes: (1) A thorough understanding of the physics of a high- 
density electron-hole gas in ZnO is lacking. (2) Around the exciton resonance at 3.31 eV, 
no pump-probe reflectivity measurements have yet been reported. (3) All experiments were 
carried out at very high carrier densities, far above tim- The ultrafast carrier and reflec- 
tivity dynamics in the exciton regime and across the crossover from the exciton regime to 
the electron-hole-plasma (EHP) regime have not yet been studied. (4) For a straightforward 
theoretical analysis of pump-probe reflectivity data one needs a homogeneous carrier density 
within the penetration depth of the reflected probe. It is the aim of the present paper to 
report in considerable detail on progress in all these four directions. This is achieved in the 
following manner. 

In Sees. HT1 and ITTT1 of this paper, we present pump-probe reflectivity data on a ZnO single 
crystal, taken at four probe wavelengths around the exciton resonance. We used 800-nm 
pump pulses to ensure a homogeneous carrier density within the penetration depth of the 
reflected probe. Excitation took place via three-photon absorption (3PA). Measurements 
were performed in a broad density range of 10 22 — 10 26 m -3 , to probe the dynamics both 
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above and below the Mott density. 

The experimental data are compared with theory. In Sec. IIVI we compute the Mott den- 
sity and the electron-hole chemical potential. Using the solutions of the statically screened 
Bethe-Salpeter equation we then compute the density-dependent optical spectra. In Sec. [V] 
by comparing our theoretical results with the experiment, conclusions will be drawn about 
the ultrafast carrier dynamics. Finally in Sec. IVI} the obtained intensity-dependent 3PA 
coefficient is tested by a Z-scan measurement. 

II. PUMP-PROBE METHOD 

For the experiment, 800-nm laser pulses from an amplified 1-kHz Tksapphire laser were 
split into a pump pulse and a probe pulse. The pump pulse was sent through a 500-Hz 
chopper wheel and a delay line, and then focused into the sample. The probe pulse was 
focused into a 4.5-mm thick sapphire crystal for self-focusing and white-light generation. The 
beam was subsequently sent through a BBO crystal for sum-frequency generation of 800- 
nm light and a selected frequency from the white-light pulse. Undesired wavelengths were 
filtered out. Then the probe was sent through a polarization rotator to obtain s-polarization 
and focused onto the center of the pump spot on the front surface of the sample. By changing 
the orientation of the BBO crystal the probe was tuned to any desired wavelength between 
360 and 440 nm, with a spectral resolution of 2 nm (FWHM). 

As sample we used an epi-polished ZnO single crystal (5 x 5 x 0.523 mm 3 ), purchased 
from MTI Corp. It is oriented in the [0001] direction, i.e., with the c-axis perpendicular to 
the plane of the wafer. The electric field of the probe was polarized perpendicularly to the 
c-axis. Experiments were performed on the Zn face of the crystal. The angle of incidence 
was 0° for the pump and 22.3° for the probe. The pump spot on the sample was 220 fim 
in diameter (FWHM), the probe 35 /im. The probe pulses reflected at the front surface 
of the sample were detected by a photodiode and a lock-in amplifier. Measurements were 
performed with (1.41 ± 0.10)-ps and (141 ± 5)-fs pump pulses (FWHM), and with 365-nm, 
370-nm, 375-nm, and 380-nm probe wavelengths, at room temperature. 
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FIG. 1: Dynamics of the reflectivity (reflection coefficient / reflectance) of a ZnO single crystal 
surface following a 1.41-ps 800-nm pump pulse, (a-d) Pump-probe measurement results at probe 
wavelengths of (a) 365 nm (3.397 eV), (b) 370 nm (3.351 eV), (c) 375 nm (3.306 eV), and (d) 380 
nm (3.263 eV). At all probe wavelengths we took the same fluence series, (e-h) Fits according to 
the Simple Model [Eq. (|3ip ]. (i-1) Fits according to the Saturation and Cooling Model [Eq. ()33|) ] . 

III. PUMP-PROBE RESULTS 

The results of the pump-probe measurements with the 1.41-ps pulses are shown in Fig. [1] 
(a-d). We observe at 365 and 370 nm that the reflectivity increases while at 375 and 380 nm 
it decreases with increasing pump fluence. For all measurements each extremum is followed 
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FIG. 2: Dynamics of the reflectivity of a ZnO single crystal surface following a 141-fs 800-nm pump 
pulse, (a-d) Pump-probe measurement results. Fluences at all probe wavelengths according to the 
legend of (f) and (j). (e-h) Fits according to the Simple Model [Eq. (|3ip ]. Note the different 
vertical scale here, (i-1) Fits according to the Saturation and Cooling Model [Eq. (}33|) ]. 

by a fast relaxation to a plateau different from the initial level. 

Figure [2] (a-d) shows the results of the measurements with 141-fs pump pulses. Here the 
carrier densities reached are much higher because of the nonlinearity of the absorption. At 
375 and 380 nm we see again a decrease of the reflectivity during 1 ps. At the highest pump 
fruence AR/R = —0.3. At 365 and 370 nm at the highest fluences the pump-probe signal 
increases, decreases, and finally increases to a plateau. 
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FIG. 3: Pump-probe reflectivity results with a long delay. Pump: 141-fs 800-nm pulses. Fluences: 
1150 J/m 2 for the 365-nm measurement and 1313 J/m 2 for the 370-nm, 375-nm, and 380-nm 
measurements. 

Pump-probe reflectivity results with long delays are shown in Fig. |3] We observe that 
the decay from the plateau to the original reflectivity takes several hundreds of picoseconds. 

When one encounters complicated pump-probe traces as the 1313 J/m 2 -result in Fig. 
El^a), one might be inclined to search for independent processes that explain the first sharp 
peak, the subsequent dip, and the rise again to the plateau. 

Actually, as we will demonstrate in this paper, these rich phenomena can be simply 
explained by increase and subsequent decrease of carrier density, at least if we correctly 
account for the highly nonlinear relation between reflectivity and carrier density. In Sec. IIVI 
we theoretically compute the optical spectra for several carrier densities and in particular 
this relation. In Sec. |V]the theoretical results of Sec. |IV]will be combined with a model for 
the carrier dynamics. The reader who is exclusively interested in the results for the optical 
spectra and the carrier dynamics can simply first consider Figs. [9] and [11] and then proceed 
to Sec. El 

IV. THEORY 

In order to elucidate how the reflectivity changes with increasing carrier density one has 
to study how the complex index of refraction changes with increasing carrier density. Li et 
al— compared their pump-probe results with a free-carrier Drude model, described in Ref. 
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10l . This model is correct for very high carrier densities where the Coulomb interaction 
between the carriers is almost completely screened. At carrier densities about 10 28 m -3 this 
is indeed the case. However, at our carrier densities of 10 22 — 10 26 m~ 3 , and at our probe 
wavelengths, the Coulomb attraction between electrons and holes does play a major role in 
the optical properties. Indeed, below the Mott density it produces the exciton resonance 
in the absorption and reflectivity spectra. As carrier density increases, screening gradually 
destroys the exciton resonance. 

Next to the vanishing of the exciton resonance due to screening, band-gap renormalization 
determines the optical properties: the band gap shrinks for increasing carrier density due 
to exchange and correlation effects. A final important effect for the optical spectra is band 
filling. We first quantitatively address these three phenomena and subsequently compute 
the density-dependent absorption and reflectivity spectra. 

ZnO has one conduction band and three valence bands, called A, B, and C. Each of these 
four bands is twofold degenerate because of the spin degree of freedom. The valence bands 
are split by the crystal field and the spin-orbit coupling: the AB splitting equals 10 meV, the 
AC splitting 44 meV.— At present it is unclear how the band-gap renormalization behaves in 
case of multiple split valence bands. Our goal is to set up a simple description of the many- 
body physics to explain our experimental data. For that purpose we make throughout 
this paper the simplification to take only the conduction band and the A valence band 
into account. Although we realize that this will affect our results quantitatively in certain 
parameter regimes of the experiment, we do not expect it to affect the physics of interest 
to us qualitatively. Transitions between the conduction band and the A valence band are 
allowed without spin-flip for the probe polarization in our experiment E _L c— ^ In this 
two-band model there is a single band gap of Eq,o = 3.372 eV.— Further, we use isotropic 
parabolic bands and quasi-equilibrium of the electron-hole gas, so that equilibrium statistical 
mechanics can be used to describe its properties. 



A. Coulomb screening and Mott density 

In this section we consider the screening of the Coulomb interaction and compute the 
Mott density um- The Mott density marks the crossover between the density regime where 
excitons exist (the exciton regime) and the density regime where they are screened away (the 
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EHP regime). It is important to pin down this value, not only for understanding the optical 
properties of ZnO, but also to know whether certain observed phenomena in ZnO, such as 
lasing, have an excitonic nature, as is frequently claimed, or not. The published values for 
the room-temperature Mott density in ZnO largely vary-^— and range from 3 x 10 23 to 
3.7 x 10 25 m- 3 . 

The physics of unscreened excitons is equal to that of hydrogen atoms. The ground-state 
binding energy is related to the Bohr radius ao by 

^° 2m r ag' ^ 
where m r = (l/m e +l/m/j) _1 is the reduced mass of the electron-hole pair. The electron mass 
in the conduction band and the hole mass in the A valence band have been experimentally 



determined to be m e = 0.28m (Ref. |22|) and mh = 0.59mo (Ref. |23|), respectively, so 
that m r = 0.19m . Here mo denotes the bare electron mass. The exciton binding energy is 
known to be 60 meV, from which it follows, in agreement with literature,— that ao = 1.83 
nm. The Bohr radius also obeys the relation 

47rh 2 e r e 

ao = — 5 , (2) 

e z m r 

from which we extract the relative dielectric constant e r = 6.56. Note that we use SI units 
throughout this paper. 

We describe the screened Coulomb interaction by the Yukawa potential 

K(X - X') = t- ^-|x-x'|/A % /gN 

VKEqE t I X — X' I 

where X s is the screening length. The derivation of the Yukawa potential needs the ap- 
proximation of static screening,- 24 that is, screening is established fast with respect to the 
Fermi frequencies of the charge carriers. This is a good approximation if h times the plasma 
frequency 

/ 6 2 Ti 

U P = \ > ( 4 ) 

is high with respect to the Fermi energies of the electrons and holes 

e F ,i = ^-( 37 ^) 2/3 > (5) 

where i stands for e (electron) or h (hole). The conditions frxj p > Ep te and ftw p > ep,h are 
both met if n < 2.8 • 10 26 m -3 . Since in our experiment the carrier density does not exceed 
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this value, we can use the Yukawa potential [Eq. (j3J)]. Note that we always consider the 
situation that the electron density is equal to the hole density, n e = nh- This is necessarily 
true for optical excitation. This density we call the carrier density n. 
In an electron plasma the screening length is given by-2 4 - 



e e r dfi e , 

where \x e is the chemical potential. We have electrons and holes however. In an electron- 
hole plasma the screening length is related to the screening lengths of the electron and hole 
plasmas according to 

a; 2 = Kl + K,l (7) 

For the Fermi-Dirac distribution at zero temperature, Eq. ([6]) reduces to the Thomas- 
Fermi screening length. In the classical high-temperature limit, the particles have a Boltz- 
mann distribution and Eq. reduces to the Debye-Hiickel screening length.- 4 We do not 
take any of these limits, but we compute A Sje and X S} h using numerically determined ideal-gas 
chemical potentials, calculated from 

where (3 = 1/(/cbT). The electron chemical potential /i e is measured from the conduction 
band edge, the hole chemical potential Hh from the valence band edge. 

Of course, the electron-hole gas is not an ideal gas. Coulomb interactions play a very 
important role. That is why we compute the screening length. The Coulomb interactions 
affect the chemical potentials via the possible presence of excitons, but in order to be able to 
compute the chemical potentials including the Coulomb interactions, one must first know the 
screening length. For the computation of the screening length the ideal gas approximation 
turns out to be sufficient, as will be demonstrated in Sec. IIVBI Note that BGR does not 
affect the screening length at a certain carrier density. The result for A s at 300 K is given 
in Fig. |H 

If one approximates the Yukawa potential by the Hulthen potential, which deviates only 
very little from the Yukawa potential, one can analytically compute the ground-state exciton 
binding energy in the case of screening as^ 4 - 1 ^ 

. [1 - a /\ s ) 2 E if X s > o , . 
K = \ (9) 
if \ s < a . 




FIG. 4: Screening length X s at 300 K versus carrier density. Note the logarithmic scales. The 
horizontal line indicates the exciton Bohr radius do- The Mott density um, i- e - the density at 
which A s = do, is equal to 1.5 x 10 24 m -3 . Excitons only exist at densities below um- 

The Mott density um is the carrier density at which X s = a,Q. At this density E s = and 
excitons can no longer exist. Our calculation shows that um = 1.5 x 10 24 m~ 3 . This value 
is lower than the values given in literature, except for the values given by Klingshirn et al— 
Apart from screening of the Coulomb attraction between electrons and holes there is 
of course also screening of the Coulomb repulsion between carriers of the same kind. The 
effects of Coulomb repulsion and the screening of it on the energy levels are captured in the 
band-gap renormalization. 

B. Chemical potential 

The electron-hole pair chemical potential with respect to the band gap \i = \i e + 
describes band filling and is an important parameter for the optical properties. A positive /i 
means population inversion. In this section we will no longer use the ideal gas model of Eq. 
OH]), but a more accurate model involving Coulomb interactions. In this interaction model 
we take into account that below the Mott density part of the electrons and holes are bound 
into excitons. 

The unbound electrons and holes obey Fermi-Dirac statistics and have the distribution 
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functions 



Ue) = . (10) 



Excitons, however, obey Bose-Einstein statistics: 

The exciton chemical potential is the electron-hole pair chemical potential /i. The energy of 
the exciton is its kinetic energy minus the binding energy. We only consider excitons in the 
ground state, so e = Skin — E s for the excitons. In this model we suppose that if an exciton's 
kinetic energy is higher than its binding energy, it immediately dissociates. 
There are four possible spin states of the exciton, 

\8,m.) e { |0,0), |1,0), |1,1) }. 

Hence the density of states of excitons is twice as large as that of electrons or holes. Only 
excitons in the states 1 , 0) and |1,0) can be created by a photon and can recombine into 
a photon without a spin-flip. In the computation of the susceptibility therefore exclusively 
the states |0, 0) and |1, 0) are to be taken into account. In the computation of the chemical 
potential, however, all four states have to be taken into account, since they all four contribute 
to the density of states at equilibrium. This gives the following relation between the exciton 
density and the exciton chemical potential 

c0 



^ ( 2K+^) )3/2 ^ deVFTWsUe)t (12) 



" = "= + ^(^'f' 2 I d£v^/i(£), (13) 



with f ex (e) given by Eq. ( flTl) . Since the electron density is equal to the hole density, the 
electron and hole chemical potentials can be calculated from the following system of two 
equations with two unknowns: 

1 2m 
2^2^ 

for i = e and i = h and with n ex given by Eq. ffl2|) and fi{s) given by Eq. fflOl) . 

The result at 300 K is displayed in Fig. [5j Also — E s and ep = EF, e + ^f.k are shown 
in this figure. For carrier densities larger than 4.8 x 10 24 m -3 the chemical potential is 
positive (population inversion) and for increasing density it approaches the Fermi energy. 
The exciton binding energy decreases with carrier density due to screening of the Coulomb 
attraction and becomes zero at the Mott density. 
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FIG. 5: The Fermi energy ef = £F,e + £f,h [Eq. ©], the electron-hole pair chemical potential //, 
and the exciton ground energy level — E s at 300 K [Eq. (JSj)], versus carrier density. 
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FIG. 6: The pair chemical potential /x at 300 K according to the ideal gas model and the interaction 
model. 

In Fig. [6] the chemical potential in the interaction model is compared with the chemical 
potential in the ideal gas model. We find that at 300 K there is little difference, confirming 
that our calculation of the screening length in the previous section is a good approximation. 
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FIG. 7: Fraction of carriers that are bound into excitons at 300 K versus carrier density. 
C. Excitonic fraction 

In principle the theory in Sec. IIV Al overestimates the screening, since Eqs. (EJ) and (J?]) 
apply to a plasma of unbound carriers. Below the Mott density a fraction of the carriers is 
bound into excitons, and screening by excitons is weaker than screening by unbound carriers. 

In order to study how large this effect on the screening is, we compute the excitonic 
fraction n ex /n from Eqs. (jTBl . The result is shown in Fig. [7J The maximum of the excitonic 
fraction is 0.14 at n = 1.9 x 10 23 m -3 . The vast majority of carriers is thus not bound at 
room temperature. 

To estimate the error made in our calculation of the screening, we make the rather extreme 
assumption that the excitons do not contribute to screening at all. Then at n — 1.9 x 10 23 
m~ 3 , where the excitonic fraction is the highest, only a density of 1.63 x 10 23 m -3 contributes 
to screening. Instead of a screening length of 5.00 nm we find X s = 5.35 nm and the excitonic 
fraction becomes 0.156. If we repeat the calculation with this new excitonic fraction, we 
get X s = 5.40 nm and an excitonic fraction of 0.159. We conclude that the errors in the 
screening length as a result of using an ideal EHP theory are at most about 8%. 

For the rest of this paper we use the screening length from ideal EHP theory, as given in 
Fig. HI and the chemical potential according to the interaction model, as given in Figs. [5] 
and HE 
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FIG. 8: Band gap renormalization: the band gap of ZnO versus carrier density. BR: Beni and 



Rice=£ with numerical results taken from Ref. 



28; VK: Vashishta and Kalia 29 ; R: Roth et al. 



30. 



BK: Banyai and Koch 2 ^; Z: Zimmermann 31 ; IA: Inagaki and Aihara 32 ; Y: Ye et a/.—; L: Lu et al— 
D. Band-gap renormalization 

Klingshirn et al— took a phenomenological approach to relate the band-gap renormaliza- 
tion to the Coulomb screening. From experiment it is known that, within errors, the exciton 
resonance does not shift with increasing carrier density. Our measurements reported in this 
paper confirm this observation. From this result one must conclude that the increase by 
screening of the exciton resonance energy is compensated by the decrease due to BGR. This 
means that the density at which the BGR equals E = 60 meV, and hence Eq = 3.312 eV, 
is the Mott density. 

This is a valid argument, but the problem remains that it is notoriously difficult to 
determine, experimentally or theoretically, the BGR in ZnO. In the literature very different 
results for the BGR have been reported. Without claiming to be complete, we show several of 
them in Fig. [HJ Where applicable we inserted Eq = 60 meV, ao = 1.828 nm, T = 300 K, and 
Eg,o = 3.372 eV.— All of these BGR results, except one, could in principle be used for the 
determination of tim- Only the BGR graph of Banyai and Koch cannot be used, because the 
use of their formula requires that one already knows the magnitude of screening. Different 
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values for the BGR lead to different values for the Mott density. Therefore we believe that 
our method for calculating um is more reliable. 

In our calculations of the optical properties we use the phenomenological BGR formula of 
Banyai and Koch,— ^ according to which BGR exactly compensates the effect of screening 
on the position of the exciton resonance: 

I E G p - E + E s if n<n M , 
Eg = < (14) 
^ E Gfi - E a /X s if n> n M - 

E s is given by Eq. 



E. Susceptibility 

In the Appendix the Bethe-Salpeter equation is derived from quantum field theory and 
solved in the statically screened ladder approximation. It is derived that the susceptibility 
of an electron-hole gas in a direct semiconductor is given by 



Xbs{u) 



k 



£xk("). (15) 



Here, d cv is the dipole moment of the valence- conduct ion band transition, L 3 is the volume 
of the (cubic) crystal, and BS stands for 'Bethe-Salpeter'. The summation over all k-states 
of the carriers is three-dimensional with a step size of 2tc/L in each direction. A factor 2 is 
included to account for the spin degree of freedom. The auxiliary function Xk(^) is implicitly 
given by the susceptibility integral equation 

Xk(w) =x2(w)(l + ^iE y ^'lXk'M). (16) 
Here, the mean- field function x°(w), depending only on the length of k, is given by 

Xl(uj) = -d cv - , ]-J^~^ h (17) 

where Sk,i = h 2 k 2 /(2rrii) are the kinetic energies of the electrons and holes, fk,% = fi(s) 
are the Fermi-Dirac distribution functions, and j(u)) is the frequency-dependent damping, 
which we discuss in more detail below. 

Finally, Vs,|k-k'| is the Yukawa potential in momentum space, i.e. the Fourier transform 
ofEq. (J3D, 

e 2 1 
K ' |k - k '' ~ k 2 + k>2-2kk'cos6 + \-^ (18) 
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where 9 is the angle between k and k'. As convention for the Fourier transform of the 
potential we use 



K(x) 



k e ik x and V t 



dxV,(x)e" 



ikx 



(19) 



We note that it is also possible to derive Eq. ( TT5|) using an equations-of-motions approach.— 
We solve Eq. f lT6|) by using a matrix inversion method, described by Haug and Koch.— 
Because of rotation symmetry we can replace K,|k-k'| by its angle-averaged 

1 



V 



s,k,k' 



d9 sin 9 K,|k-k'| 



In 



(k 2 + k' 2 + 2kk')X 2 + 1 



(20) 



Ae e r kk' ~ l(k 2 + k' 2 - 2kk')X 2 s + lJ ' 

and transform the three-dimensional summation over k' into a one- dimensional summation 
over its length k', 

s 



1 



cv ^ I k'=0,s 



(21) 



Here s is the step size of the resulting k- summation .— 
We introduce the vertex function T k (u) as 

= T k (uj)x° k (uj). 

Inserting this into Eq. (j2Ti) . we obtain the integral equation 



(22) 



T k (u) = 1 + 



2vT 2 cL 



5^ k' 2 v sAV Xfc'M rv(w) 



(23) 



k'=0,s, 



When T k is seen as a vector, this is an equation of the form 



(24) 



where 1 is the unit vector and M(u) is the matrix 



We see that 



3 



2n 2 d, 
r 3 



-k' 2 Vs,w Xl>{u)- 



(25) 



(26) 



where 1 is the unit matrix. 
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Also a background susceptibility xl produced by the lattice, including the valence elec- 
trons, should be included. This yields the following expression for the susceptibility of the 
semiconductor: 

^) = XL + £^- X>* 2 X ,H. (27) 

^ > k=0,s,... 

Here we again performed a transformation from the three-dimensional k-sum to the one- 
dimensional fc-sum and we added an upper limit. 



F. Optical spectra 

The complex refractive index n(u) = n'(u)+in"(u) is related to the complex susceptibility 
x{w) = x'i 00 ) + ix"( u ) as n{uj) = a/1 + x(w). The reflectivity (reflection coefficient) R(u>) 
of s-polarized light (the probe light in our experiment) is related to the complex index of 
refraction^ 5 (p. 422) and reads 

R=l 4aC ° S ^ (28) 

cos 2 i + 2a cos % + a/6 2 + An' 2 n" 2 

with a = (6 2 + 4n /2 n" 2 ) 1 / 4 cos[| arctan(2n'n"/6)], b = n' 2 — n" 2 — sin 2 i and i the angle of 
incidence. The absorption coefficient is given by 

a(w) = aj+ -^- L , (29) 

where c is the vacuum speed of light and ai is absorption due to crystal impurities, relatively 
very small and frequency-independent within the frequency range of our experiment. 

In order to find the density-dependent absorption and reflectivity spectra, we thus com- 
pute for each carrier density (1) xti 00 ) fr° m Eq. ffTTl) . and (2) the inverse of the matrix 

— ^i(co>), inserting appropriate values for the chemical potential, screening length, damp- 



ing and band gap renormalization. The spectra can then be obtained via Eqs. 
and fl27UH]). 

For the computation of Xk(. u )i a damping function 7(0;) is needed, representing the com 



36 a 



bined effect of carrier-phonon, carrier-impurity, and carrier-carrier scattering. In Ref. 
microscopic description of the damping due to carrier-carrier scattering is given. In order 
to stay close to the experiment, however, we here choose to work with a phenomenological 
function 7(0;). A frequency-dependence is necessary to correctly describe the Urbach tail, 
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FIG. 9: Theoretical spectra of excited ZnO at 300 K for E 1 c: (a) real part of the susceptibility, 
(b) imaginary part of the susceptibility, (c) real part of the index of refraction, (d) imaginary part 
of the index of refraction, (e) reflectivity at i = 22.3° [Eq. (j28|) ] . (f) absorption coefficient [Eq. 
(|29j) ]. The carrier densities in all graphs are, from the highest curve to the lowest curve: 5 x 10 21 
(solid), 5 x 10 23 (dash), 1 x 10 24 (dot), 2 x 10 24 (solid), 5 x 10 24 (dash dot dot), 1 x 10 25 (dash), 
2 x 10 25 (solid), 3 x 10 25 (short dot), 5 x 10 25 (dash dot), and 1 x 10 26 uT 3 (solid). In all graphs 
the gradual disappearance of the exciton resonance due to screening is visible. 

i.e. the decrease of absorption on the long-wavelength side of the exciton resonance.— In 
our calculation we take 

7 ( W ) = e {-hw+E G -E s -E a )/E a + i ( 30 ) 

This function gives the best agreement between the theoretical results for the optical spec- 
tra at the lowest carrier densities and measured optical spectra of unexcited ZnO.— £ZrM. In 
principle, at higher densities damping is stronger because of increased carrier-carrier scat- 
tering. We choose however to work with a density-independent damping in order to reduce 
the number of parameters. 

For the numerical computation we choose a step size s = 5 x 10 7 m _1 and an upper limit 
kmax — 2.5 x 10 9 m _1 . As a result our to be inverted matrix has a size of 51 x 51. We checked 
that with smaller step sizes the same results are obtained, but with a longer computation 



18 



TABLE I: Model parameters. 



Parameter Value 

d cv 4.2 x 1(T 29 Cm 

k max 2.5 x 10^ m 
XL 2.4 
fryo 50 meV 

E a 22 meV 

a/ 1.1 x 10 3 m" 1 



time. With a larger step size one obtains unphysical fluctuations in the spectra. 

The results of the computation for the complex susceptibility, the complex refractive 
index, the reflectivity at % — 22.3°, and the absorption coefficient are presented in Fig. 
[9j All low-density spectra exhibit an exciton peak. The exciton peak in the absorption 
spectrum [Fig. [21(f)] is at 3.31 eV, precisely where it should, a first evidence that our theory 
works well. The exciton resonance disappears from the spectra at densities around the Mott 
density. This is a second support for our results, or reversely, for the value of the Mott 
density that we obtained earlier. For densities exceeding 2 x 10 25 m -3 negative absorption, 
i.e. gain, appears. 

Our theory has in principle six free parameters. Table [I] shows their values. We have 
determined these values by fitting the low-density result of our theory to experimental data 
on unexcited ZnO. Our theory therefore has predictive power for higher carrier densities. 
Five parameters were determined by fitting the theoretical spectra at low density to the 
published experimental dat a 13 ' 37 " — on the linear absorption and refractive index spectra of 
unexcited ZnO, both near and far from the exciton resonance. Our value for 70 is higher 
than the room temperature damping parameters of Refs. 
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and 



40, but lower than that 



of Ref. 



38l . In Fig. (TO] the real refractive index from our model is compared with the long- 



wavelength experimental data of Refs. [37| and [l_3|. The impurity absorption coefficient cti 
was determined by a simple measurement of the transmission through our 523 /im thick 
ZnO crystal at wavelengths around 400 nm. This measurement shows that aj = 1.1 x 10 3 
m -1 , a factor ~ 10 4 smaller than the absorption coefficients near the band gap. 

In order to connect our pump-probe measurement results to theory, we calculate the 
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photon energy (eV) 

FIG. 10: Real part of the refractive index of unexcited ZnO. Our model is compared with the 
long-wavelength experimental data of Bond 37 and Jellison and Boatner ^ 
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FIG. 11: Theoretical reflectivity for i = 22.3° at 365 nm (3.397 eV), 370 nm (3.351 eV), 375 nm 
(3.306 eV), and 380 nm (3.263 eV), versus carrier density. 

reflectivity versus density at 365, 370, 375, and 380 nm. The results are shown in Fig. [TTJ 
The kink at the Mott density in some of the graphs is due to the kink in the BGR formula 
used [Eq. flU}]. 

We see at 365 and 370 nm that reflectivity rises with increasing carrier density, reaches 
a maximum, and decreases again. This behavior can be related to the pump-probe signals 
of Figs. [TJa,b) and |!__(a,b). For 1.41-ps pump pulses the carrier density remains left of this 
maximum, while for 141-fs pulses it goes beyond the top until maximum carrier density is 
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reached at the bottom of the dip in the pump-probe signal. After that point, the density 
decreases again. 

At 375 and 380 nm reflectivity monotonically decreases with increasing density. This is 
consistent with all experimental data at those wavelengths, except for a tiny peak in Fig. 

m<o. 

V. ULTRAFAST CARRIER DYNAMICS 

Our extensive experimental results on the reflectivity versus time (Figs. QJa-d), |2Ja-d), 
and [3]), combined with our theoretical results on the reflectivity versus carrier density (Fig. 
[TTT) . allow for the determination of the carrier density versus time. As we will show now, a 
concise carrier dynamics model can be found accounting for all measurements, performed at 
different probe wavelengths, both above and below the exciton resonance, at several fluences, 
both for long pump pulses (low density) and short pump pulses (high density). 

A. Buildup and decay 

Since the band gap (3.37 eV) and the exciton energy (3.31 eV) are larger than two times 
the photon energy (1.55 eV), absorption of an 800-nm pulse is a three-photon process. 3PA 
of 800-nm pulses in ZnO has been reported by He et a/.— and Dai et alM. Thanks to the 
large penetration depth, 3PA provides a homogeneous carrier density over the penetration 
depth of the reflected probe (about 50 nm). This presents clear advantage of 3PA over 
one-photon absorption. 

Following carrier buildup, we observe at all probe wavelengths, both for 1.41-ps and 141- 
fs pulses, a fast relaxation to a reflectivity level higher or lower than the initial reflectivity 
[Figs. []Ja-d) and[2]Ja-d)]. The subsequent decay to the initial level takes hundreds of picosec- 
onds (Fig. [3]), in agreement with decay times measured in time- resolved photoluminescence 
experiments.— — This slow decay is the result of radiative and nonradiative recombination 
of carriers and excitons. The remainder of this paper is devoted to extracting the carrier 
dynamics during the first 6 ps after the pump pulse. On this timescale the slow decay can 
be safely ignored. 

In line with literature we explain the fast decay by trapping of carriers into impuri- 



21 



g 



7xi(r 




, (a) 1.41 ps 


356 J/m 2 " 


6xl(f 






534 J/m 2 - 


5X10 23 






676 J/m 2 

872 J/m 2 " 


4X11T 






QQ7 T/m 2 

yy / j/m - 


3xl(f 




\ 

\ — - 




2xl(f 


//' 






lxlO 23 











2 -1 


1 2 3 


4 5 6 



2.5xlO zo 
S ZOxlO 26 



g 1.5xl(r° 

T3 



§ 1.0x10-" 

o 

5.Qxl(f 





1 (b)141fs 534 J/m 2 " 




\ 603 J/m 2 


- 


\ 854 J/m 2 " 




\ - -11 ■SO 1/m 2 

\ -1 1 »/V/ J/ 111 




\\ 1313 J/m 2 " 

V 









delay (ps) 



1 2 3 
delay (ps) 



FIG. 12: Carrier density versus time deduced from experiment using the Simple Model for (a) 
1.41-ps and (b) 141-fs 800-nm pump pulses. 

ties, such as oxygen vacancies. It is known that the density of singly ionized oxygen vacancy 
traps in a surface layer of 30-100 nm is much higher than in the interior of the crystal.— 
The fast decay therefore mainly occurs in this surface layer. This idea is supported by the 
30-nm thick surface-recombination layer found by Shalish et al— and with the observation 
of Magoulakis et al— that increased surface roughness leads to higher trapping efficiencies. 
To explain the relaxation to the plateau, we make a distinction between charge carriers near 
the surface (the surface carriers) and charge carriers in the interior of the crystal (the bulk 
carriers). Only the surface carriers are subject to fast decay. After about 2 ps all surface 
carriers have been trapped and the remaining bulk carriers produce the reflectivity plateau. 

Alternative explanations for the fast decay, like Auger recombination and stimulated 
emission, can be excluded, since at low pump fluences and low densities the fast decay is as 
prominent in the pump-probe results as at high fluences and high densities. 

Trap saturation cannot explain the reflectivity plateaus. Such an explanation requires 
a single limited trap density. From the pump-probe results of Figs. [T^a-d) and [2]^a-d), 
however, it can be found, using Fig. CD] as a gauge, that for 1.41-ps pump pulses the plateau 
is reached after a fast density decay in the order of 10 23 m~ 3 , while for 141-fs pulses it is 
reached after a decay in the order of 10 25 m -3 . 
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B. Simple Model 



The above considerations form the basis of a simple carrier dynamics model, the 'Simple 
Model'. The carrier density is taken to be homogeneous within the penetration depth of 
the reflected probe. This density n, being the sum of the surface carrier density within the 
probe penetration depth ns and the bulk carrier density within the probe penetration depth 
rig, n = n$ + rig, determines the reflectivity according to Fig. [TTJ We do not take spatial 
variations in the refractive index into account. In this model the carriers are in thermal 
equilibrium at T = 300 K. 

The mathematical expression for the Simple Model is 

dn s (t) Sa 3 I(t) 3 n s (t) 



dt 3hu t 

dn B (t) (l-S)a 3 I(t) 



(31) 



dt 3ftw 

with the initial conditions ns(— oo) = ns(— oo) = 0. Here a% is the 3PA coefficient for 
800-nm light, fiw = 1.55 eV is the pump photon energy, r is the surface carrier decay time, 
S is the fraction of the carriers within the probe penetration depth that are near the surface, 
and I(t) is the intensity of the pump pulse 

1(f) = il^MZe- t2 /( M2 ), (32) 
V27T(i 



where R(cj) = 0.105 is the reflectivity of the 800-nm pump, / is the pump fluence, and d 



measures the pulse length: 141 / \/8 In 2 = 60 fs or 600 fs. 

By fitting the Simple Model to the experimental results, we arrive at = 5 x 10~ 27 
m 3 /W 2 , r = 0.7 ps, and S = 0.8. Our trapping time of 0.7 ps is in good agreement with 
values reported in literature iLi&SS, 

Figure fT2l shows the dynamics of the carrier density n as deduced from the measurements 
using the Simple Model. Note that for 1.41 ps pump pulses the carrier density remains below 
the Mott density, while for 141-fs pulses densities in the order of 10 26 m~ 3 are reached. 

Combining Fig. [12] with Fig. [11] yields the theoretical reflectivity versus time, Fig. l(e-h) 
and Fig. 2(e-h). The agreement with the experimental results is surprisingly good, both 
with respect to the shapes of the pump-probe results as with respect to the absolute values 
of AR/R. 

Our value for is a factor 2 lower than the value of (1.0 ± 0.2) x 10~ 26 m 3 /W 2 , reported 
by He et al—, but one must note that their value was obtained for intensities / <4x 10 14 
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FIG. 13: Thermalized carrier density versus time deduced from experiment using the Saturation 
and Cooling Model for (a) 1.41-ps and (b) 141-fs 800- nm pump pulses. 

W/m 2 , while in our 1.41-ps measurements intensities of 6 x 10 14 W/m 2 and in our 141- 
fs measurements intensities of 8 x 10 15 W/m 2 are reached. At high intensities 3PA gets 
saturated and a 3 decreases with increasing intensity.— The more sophisticated 'Saturation 
and Cooling Model' described in the next section takes this effect into account. 



C. Saturation and Cooling Model 

Comparing the theoretical traces of Fig. E]^e-h) with the experimental ones of Fig. (2(a-d), 
we see that the plateau levels are in good agreement. However, the negative peaks are too 
sharp and too large. In Fig. [1] the agreement is better, but still there is some room for 
improvement: the low-fluence peaks are too small. In this section we describe a 'Saturation 
and Cooling Model' that takes 3PA saturation and carrier cooling into account, and matches 
the experimental results better than the Simple Model described in the previous section. The 
rate equations read 

dn H (t) _ a 3 [I(t)]I{t) 3 n H (t) 

dt 3hw k{n tot ) ' 

dn s (t) Sn H (t) _ n s (t) 

dt k{n tot ) r ' [M) 
dn B (t) = (1 - S)n H (t) 

dt k(n tot ) 

with initial conditions Uh{— oo) = ns(— oo) = ng(— oo) = 0. The total carrier density 
ntot = n H + n s + n B- Saturation of 3PA is described by an intensity-dependent a 3 . Like in 
the Simple Model, r = 0.7 ps, S = 0.8, and the intensity is given by Eq. fl32|) . 
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In the Saturation and Cooling Model it is recognized that carriers are created high in 
the bands and consequently do not immediately contribute to the susceptibility and the 
optical properties at optical frequencies near the band gap. Furthermore their contribution 
to screening is negligible because of their high kinetic energy. The optical properties are 
governed by the thermalized carriers only. 

In Eq. ( 133|) nn is the density of carriers high in the bands. They are created via 3PA and 
cool down to thermalized surface carriers ns and thermalized bulk carriers ns with cooling 
time k(n to t)- The reflectivity is, like in the Simple Model, determined by n = ns + ub via 
Fig. dU 

The resulting carrier density responsible for the optical response ns + % is given in 
Fig. [HI The theoretical reflectivity versus time is given in Figs. [TJ(i-l) and[2](i-l). The 
sharp peaks of the Simple Model are smoothed and the peaks have amplitudes that are 
in much better accordance with the experimental results. Note that even the complicated 
high-fluence results of Fig. E](a,b) are faithfully described by this model. 

Gu et a/.— have studied 3PA saturation of 780-nm light. They experimentally determined 
03(1) to be 



with 0:3 Gu = 1.3 x 1CT 26 m 3 /W 2 and Is = 4.4 x 10 14 W/m 2 . In our opinion this equation 
cannot be correct at high intensities, because according to it the total carrier density does not 
exceed 4 x 10 23 m -3 , while it is evident that in our experiment at least a factor 10 2 higher 
densities are reached. By fitting the Saturation and Cooling Model to our experimental 
results we have found 

«3[/] = « 3 4 + r ^ 77 |, (35) 

with a£ = 3 x KT 27 m 3 /W 2 , af = 7 x 1(T 27 m 3 /W 2 , and I s = 4.4 x 10 14 W/m 2 . The value 
obtained for Is is the same as found by Gu et al. In the limit of low intensities a% is equal 
to the result of He et al.—. 

The cooling time k(n tot ) is density-dependent.— Carrier cooling times in the range of 30 
fs to 1.75 ps have been reported in literature.— ^ 4 -~— By fitting the Saturation and Cooling 
Model to our experimental results we found 

k(n m ) = k + rn t ot, (36) 
with k Q = 0.2 ps and r = 4 x 10~ 39 m 3 s. 
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z-position of the crystal (mm) 



FIG. 14: Z-scan. The diamonds indicate the measured 800-nm transmission through the 523 fim 
thick crystal versus the position of the crystal with respect to the laser focus. The line is the 
transmission calculated using the intensity-dependent three-photon absorption coefficient found in 
Sec. E22 i.e. Eq. ([35]) with a^ = 3x 1(T 27 m 3 /W 2 , af = 7 x 1CT 27 m 3 /W 2 , and I s = 4.4 x 10 14 
W/m 2 . 

It is interesting to compare our results with the 266-nm pump-THz probe results of 
Hendry et al— They also found a fast initial decay of the carrier density, followed by a 
plateau. In their analysis, the decay was attributed to Auger recombination. 

VI. Z-SCAN MEASUREMENT 

In order to test the 800-nm 3PA coefficient obtained in Sec. IV C] an open-aperture Z-scan 
measurement was performed. The open-aperture Z-scan has been reported for the first time 
by Sheik-Bahae et a/.— as a sensitive technique to measure non-linear absorption coefficients. 

For this measurement, the ZnO crystal was moved along the 800-nm beam through the 
focus (i.e. in the z-direction, hence the name 'Z-scan'). The absolute transmission was 
measured as a function of the position of the crystal. The results are shown in Fig. [TH For 
this measurement, the pulse duration was 135 ± 5 fs (FWHM), the maximum fluence 355 
J/m 2 , and the maximum intensity 2.2 x 10 15 W/m 2 . The Rayleigh range was measured to 
be 3.6 mm. 

The line in Fig. [TH is the transmission through the ZnO crystal, calculated from the 
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intensity-dependent 3PA coefficient 0:3 [J] stated above [Eq. (155]) ]. and also taking into 
account the 800-nm reflectivity at the front- and backside of the crystal. Evidently, the 
agreement between the measured and calculated transmission is excellent. This result con- 
firms the value for the 3PA coefficient found in Sec. IV CI Thus, it is also further evidence for 
the reliability of the theoretical spectra shown in Figs. M and [TTJ and the carrier dynamics 
shown in Fig. [131 

The transmission data of Fig. [TH also confirm that the penetration depth of the 800-nm 
pump pulses is very long compared to the wavelength of the probe. Therefore the 800-nm 
pump indeed creates a homogeneous carrier density within the penetration depth of the 
reflected probe, in contrast to a pump with a photon energy above the band gap. 

VII. CONCLUSIONS 

The ultrafast screening and carrier dynamics in ZnO, including the crossover between 
the exciton regime and the electron-hole plasma regime, have been studied experimentally 
and theoretically. Pump-probe reflectivity measurements, taken at four probe wavelengths 
near the exciton resonance and in a broad range of carrier densities (10 22 — 10 26 m -3 ), show 
rapid (~ 1 ps) and strong (up to 30%) changes in the reflectivity. These effects result from 
the disappearance of the exciton resonance due to screening. Other processes affecting the 
optical properties at high densities are band-gap renormalization and band filling. 

Our calculations show that the Mott density in ZnO is 1.5 x 10 24 m~ 3 at 300 K. This 
means that phenomena occurring at higher carrier densities in ZnO cannot be related to 
excitons. This includes lasing in ZnO nanowires and other ZnO nanostructures: if the 
density is higher than 1.5 x 10 24 m -3 , lasing must be electron-hole plasma lasing.— 

To calculate the optical spectra of highly excited ZnO, we used quantum field theory of 
a quasi-equilibrium system of electrons and holes that interact via the screened Coulomb 
potential. We computed the density-dependent spectra of the complex susceptibility, the 
complex refractive index, the reflectivity, and the absorption coefficient. By comparing 
the theoretical reflectivity spectra with the pump-probe data, we obtained a description 
of the carrier dynamics, consisting of 3PA with an intensity-dependent 3PA coefficient, a 
density-dependent carrier cooling time, and a distinction between surface carriers having a 
decay time of 0.7 ps and bulk carriers having a decay time of hundreds of picoseconds. The 
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agreement between the theoretical reflectivity based on this model and the experimental 
results is excellent. Finally, the results of an open-aperture Z-scan confirm the obtained 
3PA coefficient. 

These results provide strong evidence that this many-body theory well describes screening 
and band filling in ZnO at high carrier densities and that the band-gap renormalization is 
approximately given by Eq. (TH|) . The crossover between the exciton regime and the electron- 
hole plasma regime, as well as the optical spectra, are faithfully described by this theory. 
Our results for the density-dependent optical spectra and the ultrafast carrier dynamics are 
of general importance for ZnO research. 
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Appendix 

For the computation of the optical properties of a semiconductor, an expression for the 
susceptibility is needed. In this Appendix we derive the susceptibility of an electron-hole gas 
in a direct-band-gap semiconductor from many-body quantum field theory within random- 
phase approximation (RPA). In particular, we derive Eq. ( fl5|) from the statically screened 
Bethe-Salpeter ladder equation. 

1. Polarization and susceptibility 

Consider a direct-band-gap semiconductor crystal, subject to an oscillating external elec- 
tromagnetic field with angular frequency u, 

£(t)=£ e-^. (A.l) 

In the pump-probe experiment of this paper £(t) is the field of the probe laser pulse. The 
electromagnetic field incites a polarization response of the electron-hole gas present in the 
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semiconductor: The polarization oscillates at the same frequency, 



P{t) = P e- lwt . (A.2) 

As long as the electric field is not extremely strong, there is a linear relation between Sq and 
Pa, 

P = x(u)e £ L 3 . (A.3) 

The polarization response is thus described by the complex susceptibility xi 00 ), where we ex- 
plicitly indicate the dependence of the susceptibility on the frequency of the electromagnetic 
field. Our probe laser pulses were weak enough to be in this regime of linear response. 

When after the pump pulse the charge carriers have relaxed to Fermi-Dirac distribu- 
tions at a certain temperature, equilibrium statistical mechanics can be used to describe its 
properties. The expectation value for Pq is given by the relation 

where Z gr is the grand canonical partition function. Computing this Z gr , we use the many- 
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body quantum field theory described in chapter 12 of Ref. 



Action 



The grand canonical partition function describing the response of the electron and hole 
fields in the direct-band-gap semiconductor to the external field S is given by the functional 
integral 

Z gr (£) = [ d[0*]d[0]e- 5 ^' £ *<^ (A.5) 
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where stands for <p e f, e |, 0^, 0/^, the electron and hole fields with spin up and spin down, 
respectively, and S is the action, given by 



Here a stands for "| or 4, i again stands for e or h, r = it is imaginary time, and the Yukawa 
potential V s is given by Eq. 

The first term describes the energy of the quasi-free electrons and holes (quasi-free because 
the interactions between the carriers renormalize the band gap). The second term describes 
the attractive interaction between electrons and holes with opposite spin. The repulsive 
interactions and the attraction between electrons and holes with parallel spin are contained 
in the BGR and therefore not represented here. The third term describes the annihilation of 
electron- hole pairs by the electric field (stimulated emission). The fourth term describes the 
creation of electron- hole pairs by the electric field (absorption). We consider only transitions 
without spin-flip, so that the hole of the electron-hole pair always has a spin opposite to the 
electron spin. 

In the rest of this Appendix, we write fg(x) for the product of functions f(x)g(x). In 
addition, in order to simplify the theory, we replace the interaction potential V s (x — x') by 
a point interaction 



In the end, we again replace the point interaction by the screened Coulomb potential. 



S[cf>*, </>,£*,£} 




(A.6) 



x Vs(x - x!)(j) h ,- a (x', r)0 e)Q (x, r) 




K(x-x')->_M(x-x'). 
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3. Hubbard-Stratonovich transformation 



Now we introduce two fields, A-f(x, r) and A^(x, r), of which the averages are given by 

(A t (x,r)> =V (M er (x,T)), (A.7) 

(A^(x,t)> =V (<l> ht Mx,T)). (A.8) 

With these two fields we perform a Hubbard-Stratonovich transformation. The following 
two identities are inserted into the integrand of the partition function, Eq. ( 1A.5j) : 



l = exp[Tr[ln(-K ) - 1 /^)]] / d[A^]d[A 



r 1 P P f 

x exp {- J drydx[A;(x,r) - C^Jx, r) V ] ( A ' 9 ) 
x V r ~ 1 [A a (x, r) - V (j) h _ a eja (x, r)] 



for a =|" and a =\,. This cancels the fourth-order term in the action, so that the 
fermionic integrals become Gaussian. In the following computation we absorb the factor 
exp[Tr[ln(— V^ _1 /fr)]] into the integration measure. 

The resulting action, only quadratically depending on the fermionic fields, is 

rh$ r I A I 2 -i- I A I 2 
S[A*,A, ( p*^,£*,£] = - j drydx 1 t! ±± + ' 

rW r 

- hS j drdr7dxdx'0* Q (x, r)G Q /(x,^ x'rO^x', r) 

- ^ J dry dx[d cv £*(r) - A* (x, r)](j) h _ Q ejQ (x, r) 



(A.10) 



~Y1 / dr dx l d cv£(r) - A a (x,r)]0* a ^_ Q (x,r), 

where A stands for A^, A^ and where (x, r; x', r') and Cr^ (x, r; x 7 , r') are the inverse 
noninteracting Green's functions, given by 

G^x,r, x y) = -i{ft|- - ^! - ^}5(x - x')5(r - /). (A. 11) 
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4. Integration over the fermion fields 



The integrand of the action can be written in the form of a matrix multiplication as 



S[A*, A, <f)*,<f), £*,£} = -J drjd 



hp 

h I drdr' / dxdx' 
/ 



V <M x - f) J 



|A t | 2 + |A 4 
V 



/0et(x',r') 



(A.12) 



V^ t (x',T')/ 



Here, G 1 is the inverse Green's function matrix, which can be expressed in a noninter- 
acting part and a self-energy part 



G-^x, r; x', t') = Go A (x, r; x', r') - E(x, r; x', r'), 



where the noninteracting part G 1 (x, r; x', r') is given by the matrix 

/G^for^r') 

-G , 1 (x',r';x,r)0 
C^for^r') 

0-G„ A 1 (x' J T , ;x J r)y 



(A.13) 



V 



(A.14) 



and the selfenergy is 



S(x,r;x',r') = ^(x-x')5(r-r / ) 



/ 



x 



A t (x,r) 



A|(x, r) 



V 



o 
o 








A*(x,r 






A 4 (x,r 




(A.15) 



where A a (x, r) = A Q (x, r) — d cv £{r). 

We can now perform the integration over the fermion fields, using the well-known results 
for Gaussian integrals.— (p. 28) The result is the effective action 



S cS [A:A,S:S} = -j dr^dx 



I At 



IA, 



V 



nitpnC-G- 1 )], 



(A.16) 
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related to the partition function as 

Z gr =Jd[A*]d[A}e- scB ^' A ^^ h . (A.17) 

The trace is to be taken over space, imaginary time and over the 4x4 matrix structure of 
the Green's function (Nambu space). 

5. Power expansion in A(x, r) 

Now we expand the effective action into powers of A(x, r). We write 

G-^Gq^I-GoS) (A.18) 
and make a Taylor expansion of the logarithm, 

CO 

Tr[ln(-G- 1 )]=Tr[ln(-G 1 )] - £ -Tr[(G S) m ]. (A.19) 

m=l 

The first-order term (m — 1) is 0: There is no contribution in the effective action that is 
first-order in A(x, r). For the second order term (m = 2) we find 

h l i h K 



2 



1 f p f 
Tr[(G S) 2 ] = -j-Ysj dTdT 'J dxdx ' 



(A.20) 



x G 0e (x, r; x', r')A Q (x', r)G oh (x, r; x', r')A;(x, r). 

Now we make the approximation to ignore all terms in the effective action with order 
> 2. Within this approximation the effective action can be written in matrix multiplication 
form as 

S cS [A*,A, £*£} = -ftTrM-Go 1 )] -J drdr' J dxdx' 
±55+\G 0e G oh 




+ \G()eGoh 



\Go e Gohd C v£(j r ) 

\GoeGohd C v£(j') 

— \ Goe Goh d cv £ ( T ) 

—\Go e Gohdcv£(T~) I \Aj 

2 „ 1 

+ ■^Go e Gohd cu £*{T)d cu £{T ) >, 
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(A.21) 



where (A|,A£) is a shorthand for (Aj(x, r), A|(x, r)), (A t , A 4 ) for (A t (x', r'), A ; (x', r')), 
W for 5(x - x')<5(r - r'), and GoeGoh for G 0e (x, r; x', r / )G h(x, r; x', r'). 

6. Bethe-Salpeter ladder equation 

In the normal phase the expectation value for the A-fields is zero. It only becomes 
nonzero if the temperature decreases below a certain critical temperature, far below room 
temperature. If that happens, condensation of electron-hole Cooper pairs occurs and the 
system becomes a superfluid.— Here we describe the optical properties in the normal state. 
Since the effective action up to this order is quadratic in the A-fields, the integration over 
them can easily be performed:— (p. 26) 

Z gr (S) = e - s[£ *' £]/ \ (A.22) 

where S[£*,£] equals, including the factor absorbed into the integration measure, 

S[£*,£) = S BS [£,£*} -2hTr[ln(-V - 1 /h)) 

( ±66 + ±G 0e G oh (A (A.23) 
-/m-fh^-Go 1 )] +Tr[ln(- I h )]. 

\0 ±55+lG 0e G 0h J 

Note that the last term was obtained by integration over boson fields, while the third term 
resulted from integration over fermion fields. This explains the opposite signs. 

The term Sbs[£*, £] is the Bethe-Salpeter light-matter action. Being the sum of all action 
terms containing £(t), it determines the optical properties. From now on, we concentrate 
purely on this term. The action S^g [£*,£] is given by 

S BS [£*,£] = S MF [£*,£] + S CA [£\£\. (A.24) 

The term Smf[£*,£] here is the mean-field action 

S MF [£*, £} = -J o drdr'Jd^d 2 cv £*(r)-G 0e G oh £(r'). (A.25) 

The mean- field action is the part of the action without Vq, so without the attractive interac- 
tion between electrons and holes with opposite spin. It is like a free-particle term; we call it 
a 'quasi-free-particle term', since the BGR includes the other Coulomb interactions. Figure 
IT5Ta) shows the Feynman diagram of the mean-field action. 
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FIG. 15: Feynman diagrams of the Bethe-Salpeter light-matter action: (a) Mean field, (b) Electron- 
hole Coulomb attraction. 




= V\/\/V + G o, 




K/WVT 



V/\/\/M + 
WW 



- \zwv + 




Vwv* 



FIG. 16: Bethe-Salpeter ladder equation: Expansion of the many-body T matrix into ladder 
diagrams. 



The Coulomb attraction term, Sca[£*, £], describes the effect of the attractive interaction 
between electrons and holes with opposite spin on the optical properties. The integration 
over the A-fields shows that it equals 



drdr'dr"dr"' / dxdx'dx'dx'" 



Sca[£*,£] 



x d cv S*(r)-G 0e (x, r; x", r")G oh (x 7 r; x", r") 

x(^+kGo^ 1 (xV";x"'r"') 
v Q n 

x \G, e {^\ r'"; x', T')G oh (x"', t">; x', r')d cv £{r') 
a 



(A.26) 



Eqs. ( 1A.241lAT26l) constitute the Bethe-Salpeter ladder equation. 

Figure IToT b) shows the Feynman diagram of the Coulomb attraction term. The square 
denotes the many-body T matrix. It consists of the ladder sum of diagrams shown in Fig. 
TBI The interaction V s here is the statically screened Coulomb potential (Yukawa potential). 
By performing a little algebra with the diagrams, as shown in Figs. [16] and [171 we obtain 
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( 1 - ^"j^^l 00 " ) = => 




FIG. 17: Derivation of Sca\£* ■> £] [Eq. ()A.26p ] via Feynman diagrams. 



a result for the many-body T matrix that can be recognized in Eq. ( 1A.26I) . where a point 
interaction replaces the Yukawa potential. We see here that the Bethe-Salpeter light-matter 
action can be derived in two ways: analytically in the way explained above, and alternatively 
via Feynman diagrams. 



7. Fourier transformation 

In order to find the susceptibility as a function of frequency the light-matter action has to 
be transformed from coordinate space to momentum space and from imaginary time space 
to imaginary frequency space. We start by transforming the Green's functions. 

The inverse Green's functions are defined by Eq. (lA.llj) . By using the definition of the 



inverse of a matrix in coordinate and imaginary-time space, 



6{x - x>(r - t') 

f hf3 r (A. 27) 

= J dr'J dx'M a , Q ,(x, r; x', t')M~^ , r'; x", r"), 



we find for the Green's functions 

1 °° _fc ik-(x-x') p -to„(T-T') 

G ,(x,r; x', t') = -^CE * / > (A.28) 

hpL d ^^ -ihuj n + e ki - fii 

k n=— oo ' 

where £k,i = £fc,i = h 2 k 2 /(2nii), and u n are the fermionic Matsubara frequencies u n = 
ir{2n + l)/(hp). 

For Fourier transforming the screened Coulomb potential we use Eqs. (fT9l . but for other 
functions we adopt the conventions 



oo 

,(v, ^?l /(k * ,t4Y ~ ; 



1 r n P r 



(A.29) 



Further, 

^-Vk.(x-x') = L 3 5(x - x'), and J ^e^" 1 *')* = 5 kk ,. (A.30) 
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The Fourier transform of Eq. (1A.28j) is given by 

-h 

G i(k,iu n ]k',iuj n >) = 5k.k'5„,„'— ■ (A. 31) 

-ihoj n + e kii - fa 

or, in a shorter notation, 

G ,(k, iu n ) = — — . (A.32) 

-ihw n + e kti - Hi 

Being a photon field, the oscillating electromagnetic field is expressed in quantum field 
theory in bosonic Matsubara frequencies u)f — 2ixf/ (hj3), and in imaginary time r = it, 



£{t) = - 7 =£ f e^f T , and £*{t) = -j=£}e w f\ (A.33) 



=£ f e^f T , and £*(r) = -L. 

where Sf = y/Wfi £q. The relation between the Matsubara frequency and the real frequency 
is given by the Wick rotation ihujf — >■ flu — Eg + He + Hh- Because we work here with 
complex electric fields instead of real fields, an extra factor 1/2 has to be included into the 
terms quadratic in £. 



8. Mean-field susceptibility 



Inserting Eqs. ( 1A.28|) . ( 1A.32j) . and (1A.33j) into Eq. ( 1A.25j) . and elaborating the integrals 



over coordinate space and imaginary time, we obtain for the mean-field action 
S M f[£*,£] = --^^2^2d 2 cv £* f G 0e (k,iLj n ) 



' k,k' n,ri 



X 



Cofc(k , iw n /)fj/5_k,k'5k,-k'5 

= -^2pd 2 cv \£ f \ 2 ^ ^2 Gr ° e ( k ' iu n )G oh (-k, iuf- n -i) (A.34) 

k n 

= \p 12 V — 

h2(3 cv\ /I Z^Z^_ ihw + 

k n ' 

-h 

where we used Eq. (1A.30|) . The Kronecker deltas describe the conservation of momentum 
and energy. Figure [T8l(a) shows the corresponding Feynman diagram. The momentum of a 
created electron is equal and opposite to the momentum of the created hole. The relatively 
very small momentum of the absorbed and created photons is neglected in this derivation. 
The third and fourth Kronecker deltas imply / = n + n' + 1 = /', meaning that if an 
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FIG. 18: Feynman diagrams of the Bethe-Salpeter light-matter action in momentum and Matsub- 
ara frequency space: (a) Mean field, (b) Simplest diagram of the electron-hole Coulomb attraction. 

electron-hole pair created by a photon recombines again, the emitted photon has the same 
frequency as the absorbed photon. 

We split the fraction and perform the sum over n by contour integration, 62 (p. 140) 

Smf[S\ £] = lim V — 



x 



+ - (A.35) 



-ihLO n + E^e - fl e -itih}f- n -l + £ k;h ~ Hh 

= d\ v ml E - k — 1 ~ /k,e " /m , 

^ ifrhjf - £ ki£ - £ kj/l + He + Hh 

where f^ e = f^ e and /k,/i = fk,h are the Fermi-Dirac electron and hole distribution functions. 
Wick rotating to real frequencies and including a frequency- dependent damping factor, 

ihcjf —> H(oj + — Eq + /i e + Hh, (A. 36) 

we obtain for the mean-field action 

S MF [S*,£] = -d^ElY.xl^), (A.37) 

k 

with 

x&O = x» = -d c l -f^-f^ (A . 38) 

ti{u + ij{u))-ek,e-£k,h-EG 
This mean-field function can be used to compute the optical properties within mean-field 
theory, i.e., by only taking into account the mean- field action. The mean-field contribution 
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to the partition function is 

Zf/{£) = eS^M* = exp[4,/3^ E Xk(«)] (A.39) 

k 

which, using Eq. ( 1A.4j) . gives a mean- field polarization of 

< Pq >mf= 2d cv £o Xk(w), (A.40) 

k 

so, according to Eq. (1A.3I) . the mean- field susceptibility equals 

X^M = ^j£x&0- (A.41) 

9. RPA susceptibility 

In order to obtain the susceptibility from the full Bethe-Salpeter ladder equation (i.e. 
within RPA), we need to transform also the Coulomb attraction term [Eq. (1A.26I) ] into 
momentum and Matsubara frequency space. 

First the factor {^SS + \GoeGoh) 1 has to be expanded. From Eq. (IA.27I) we see that 

(^SS+lG 0e G oh )-\K',T'^",r") 
V n 

= V 5{x' - x>(t' - t") 

r - Co e Go/i(x , r ; x ,r ) 

n 

v 3 r h/3 r 

+^ J dr'ydx'"Go e G ,(x',r';x w ,r w )G'oeGo/ ^ (x'V" / ;xV /, ) (A.42) 

T/4 [hp f 

' drV" / dx w dx w/ G 0e G r o? l (x / , r'; x'", r w ) 



xGo e Go,(x'V w ;x ,/ V ,/,/ )G'oeGo /l (x' w ,r' w ;xV /, ) 
+ . . . 

In this expansion the ladder sum of Fig. [16] can be recognized. 

By inserting Eq. flAT42|) into Eq. flAT26|) . and by using again Eqs. f lA~28]) . f lA~32]) . flA~33|) . 
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and (1A.30j) . we obtain for the Coulomb attraction term 

S C a[£*,£] = ^p d cv\^f\ 2 

k,l n,m 

x G 0e (l,ici; ni )Go&(-l,iw/_ m -i) 
~ ¥¥T 6 E G r oe(k,ia; n )Go /l (-k,ia; / _„_i) ^ A ' 43 ^ 

k,l.m n,m,p 

x G? 0e (l,iw m )Croh(-l > itt;/- TO -i) 
x GQ e {m,iu p )G oh {-m,iuj f - p -i) 

+ ...}. 

Figure fT8T b) shows the Feynman diagram corresponding to the first of these terms. Again 
the principles of conservation of momentum and energy follow from the theory. Every next 
term is equal to the previous one, multiplied by 

-V 



IWf—n—l)- 



k n 

Apart from a prefactor, this factor is equal to the mean-field action [see Eq. (1A.34I) ]. It can 
therefore be written as 

k 

It directly follows that the Bethe-Salpeter light-matter action, including the whole Bethe- 
Salpeter ladder, equals 

S BS [E*,£] = -d cv h(3£^J2x^), (A.44) 

k 

where 

Xk(w) = X» (l - £|l E M")) ■ (A-45) 
Finally we have to replace the point interaction by the screened Coulomb potential, 

-Vb5(x-x') -»■ K(x-x'). 

The expression for the susceptibility demands a potential in momentum space. From Eqs. 
(|T9|) and (1A.30j) it follows that in momentum space Vq has to be replaced by — V s ,ik-k'|, so 



Xk(w) = x2(w) (l + ^3 E K,|k-k']Xk'(a;)) . (A.46) 
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Using Eqs. (1A.22j) . (1A.4j) . and ( 1A.3j) . one then easily finds that the susceptibility of the 



electron-hole gas is given by 



k 



£xk(w). (A.47) 
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Note that in Ref. 



24l the factor s/(27r) 3 from the new step size, and the factor Airk from the 



integration over the angles, are not explicitly shown. 
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